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Abstract 

Laboratory-based optical analogs of astronomical objects such as black 
holes rely on the creation of light with an extremely low or even vanishing 
group velocity (slow light). These brief notes represent a pedagogical attempt 
towards elucidating this extraordinary form of light. This paper is a contribu- 
tion to the book Artificial Black Holes edited by Mario Novello, Matt Visser 
and Grigori Volovik. The paper is intended as a primer, an introduction to 
the subject for non-experts, not as a detailed literature review. 
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1 Motivation 



Creating a black hole with humble human resources is certainly a fantastic idea, 
yet perhaps not entirely lunatic. Recent experimental progress in quantum optics, 
in particular the generation of slow light [|l], |^, || has opened a route towards the 
formation of optical event horizons for light or towards other less dramatic but 
equally interesting phenomena. Most of the ideas have been discussed in a series 
of papers 0] and review articles The underlying atomic and optical physics is 
perhaps less familiar to the audience of this book. Therefore, it might be worthwhile 
to develop some of the physics behind slow light from basic principles that are 
hopefully known to most readers. 



2 Light-matter interaction 

Consider an atomic medium capable of interacting with light, for example a cell filled 
with Rubidium vapor or one of the alkali Bose-Einstein condensates. The atoms that 
constitute the medium are electrically neutral but polarizable by light. Since atoms 
are usually much smaller than an optical wavelength, light experiences the atoms 
as dipoles (as the lowest order in the multipole expansion of an electrically neutral 
charge distribution 0]). In the dipole approximation, the energy density of the 
light-matter interaction is given by the negative scalar product — P ■ E of the matter 
polarization P (dipole-moment density) and the electric field E 0. Therefore, the 
Lagrangian of light and atomic matter reads in SI units 

^ = I _ ^2^2^ +p.E + ^A. (1) 

The first term describes the free electromagnetic field with c denoting the speed of 
light in vacuum. The last term of the Lagrangian characterizes the internal dynamics 
of the atoms. We represent the electric field E and the magnetic field B in terms of 
the vector potential A in Coulomb gauge 

E = -dtA, B = VxA, V-A = 0. (2) 

Let us assume, for simplicity, that both the electromagnetic field and the medium 
are uniform in two spatial directions in Cartesian coordinates, but may vary in the 
third direction z. Furthermore, we consider light with fixed polarization so that we 
can concentrate on one component A of the vector potential A. The Lagrangian 
simplifies to 

^ = I {{dtAf - c\d,Af) -PdtA + ^A. (3) 
We obtain from the Euler-Lagrange equation the wave equation 

{d^,-c-dt)E = -^dlP. (4) 

Atoms have a well-defined level structure such that light, oscillating 10^^ times per 
second, must match the atomic transition frequencies, because otherwise the effect 
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of E in — P ■ E is rapidly washed out. Given a certain frequency range of light, the 
optical field thus interacts with a selected few of the atomic levels, which greatly 
simplifies matters. 

Let us describe the atoms quantum-mechanically, while regarding the electro- 
magnetic field as classical. The simplest relevant atomic system involves just two 
levels, say the ground state | a ) and the excited state \b). The Hamiltonian of the 
two-level atom is simply 

H = huJat\b){b\-^{A + A^)E, (5) 

A = \a){b\, (6) 

where ujab denotes the atomic transition frequency and Kab corresponds to the atomic 
dipole moment, a real number for transitions between bound states. In the Heisen- 
berg picture the transition operator A oscillates with positive frequencies near Uab- 
Therefore, A couples entirely to the negative-frequency component E^~^ of the elec- 
tric field, while couples to the positive- frequency component E^~^\ We arrive at 
the effective Hamiltonian 

H = hu,,\b){b\-'^{AE(-^ + A^E(+^) . (7) 

To describe the quantum state of the atom, we employ a density matrix p that 
characterizes a statistical ensemble of pure states | tpa ) with probabilities pa, 

P = IIPa|V^a)(^a| • (8) 
a 

Probabilities are non-negative and sum up to unity. Consequently, the density ma- 
trix has non-negative eigenvalues and is normalized as trp = 1. In the Schrodinger 
picture, the density matrix evolves while the operators are invariant in time. Ac- 
cording to Lindblad's theorem |^, the evolution of a normalized and non- negative 
density matrix is governed by the master equation 

^ = 1 [p, ^] - ^ 7. [A] Ap-2ApAl + p A] Ai) . (9) 

The Lindblad operators Ai describe dissipative processes occurring at the rates 7;, 
for example the spontaneous emission from the excited state to the ground state. As 
a result of the light-matter interaction, a medium of riA atoms per volume generates 
a matter polarization of 

P = nA'^tT{p[A + A^)}=nAtiabRe{a\p\b). (10) 

In this way the response of the atoms to the electric field modifies the propagation 
of light given by the wave equation (1). 

Consider a medium at rest in a regime of linear response. Here the medium 
integrates the local history of the electric field via the susceptibility 

/ + 00 
x{t-t')E{t')dt'. (11) 
-oo 
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Causality implies that P must not depend on the future of the E field, which restricts 
the integral (|ll|) to the time interval (—00, 0] by requiring 



X{t) = for t<0. (12) 
Consider the Fourier-transformed (spectral) susceptibility 

Xiuj) = i- xit) e'-' dt = i- xit) e^-* dt . (13) 
Ztt j~oo ztt Jo 

Because x(t) is real, we get 

x{-uj) = r{^)- (14) 

Let us regard x(^) ^ function of complex frequency uj. When the imaginary part 
of u is positive, x{^) cannot have singularities, because here the Fourier integral con- 
tains a factor exp(— [Ima;]t) that enforces convergence. Therefore, x{^) is analytic 
on the upper half plane. Causality thus imphes analyticity ^ . For a non- dispersive 
medium x{^) is constant over the relevant frequency range and x(t) is reduced to a 
delta function, describing an instant response of the medium. In dispersive media, 
x{^) varies and the poles of x{^) on the lower half plane correspond to atomic res- 
onances. At the real co axis the real part of x describes the dispersive properties of 
the medium, whereas the imaginary part accounts for dissipation. Given an analytic 
function x{w) on the upper half plane that decays sufficiently fast when u — > 00, 
the real and imaginary parts of x{^) the real uj axis are related to each other by 
Hilbert transformations (Kramers-Kronig relations 0). The imaginary part of 



x(a;) is thus uniquely determined by the real part and vice versa. 

Dispersion influences the group velocity of light pulses. Suppose that the medium 
properties do not vary significantly within the scale of an optical wavelength. In this 
case we can characterize completely the propagation of light pulses by the dispersion 
relation between the wave number k and the frequency u, 

fc' = ^ [1 + Xi^,z)] , (15) 

derived from the wave equation by Fourier transformation with respect to space and 
time. A light pulse propagates like a particle with Hamiltonian u and momentum 
k, subject to Hamilton's equations 

dk du) dz doj 

dt"~9J' ~dt^M' ^ ^ 

The group velocity Vg is the velocity dzjdt of the fictitious light particle. 



2n duj 

Slow light 0] involves an extremely dispersive medium where dx/d^ is large and, 
consequently, where the group velocity is very low (a few meters per second, or 
lower) . 



4 



3 Ordinary media 



Before we embark on discussing extremely dispersive media, let us consider an or- 
dinary medium composed of two-level atoms at rest. Assume that the atoms are 
identical, with equal transition frequencies Uab, and that they are affected by dissi- 
pative relaxation processes. The dissipation transfers excitations from the excited 
to the ground states, described by the Lindblad operator (j^). Assume that the 
transition rate 7 dominates the time scale of the light-atom interaction. In this 
regime, relaxation forces the atomic dipoles to follow the fields. In the case of linear 
response we can decompose the electric field into Fourier components. To analyse 
the response, it is sufficient to study the reaction of one of the atoms to a single 
monochromatic wave with frequency 00. We characterize the field-strength of the 
wave in terms of the Rabi frequency Q with 

Qe-''^' = ^EM. (18) 

In the absence of relaxation, an atom would oscillate between the ground and the 
excited state with frequency Q (Rabi fiopping). On the other hand, relaxation leads 
to a stationary state. To describe the stationary regime we use an appropriate 
interaction picture. 

In an interaction picture, indicated by tildes over operators, the dynamics with 
respect to a partial Hamiltonian Hq is separated from the total evolution of the 
density matrix, 

p = U^pUo, f/o=exp(^-^i^ot) ■ (19) 

To derive the evolution equation, we differentiate p with respect to t and apply the 
master equation (^. Suppose that the commutator between Ai and Hq is propor- 
tional to Ai, 

[Ai,Ho]=nujiAi. (20) 

In this case UqAiUq gives Ai exp{—iujit) , and thus the dissipative part of the mas- 
ter equation (|^) remains the same in the interaction picture. The Hamiltonian is 
transformed according to 

H = U^oHUo-Ho, (21) 

such that we obtain 

^ = ^ [p. ^] - {^l P - P^] + p Al Ai) . (22) 

Returning to the two-level atom, we use the interaction picture with respect to 

Ho = huj\b){b\ (23) 

that preserves the dissipative dynamics and leads to the time-independent Hamil- 
tonian 

H = h{ujab-uj) \ b){b\-UAn* + A^n) . (24) 
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The atom reaches a stationary state when the relaxation balances the optical tran- 
sition, 

^ [p, if] = 7 (it i p _ 2i p it + p . (25) 

Given the normalization trp = 1, we can easily solve the linear equation (p5D for the 
density-matrix components. Assuming linear response, we linearize the solution p 
in the Rabi frequency Q and obtain 
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2{UJ- UJab - il) 





(26) 



2{UJ - UJab + il) 

The positive-frequency component of the matter polarization (|l^) is therefore 



p{+) = ^^^^^(^l)\(j^plJ^^\a) = '-^ KahTab^e 



-iujt 



(27) 



where we have introduced 



Tab 



1 



1 



UJ -UJab + il 



Ua 



[28) 



For monochromatic light, p(+) gives simply eqx^^'^^- Considering the property (H) 
of the spectral susceptibihty we obtain the Lorentzian 



X 



nA f^ab 
4 Eoh 



[Tab + Ua) 



(29) 



Transforming to a real-time susceptibility x(t) we see easily that x{t) responds 
within the relaxation time 7. In accordance with causality, the Fourier-transformed 
susceptibility x is analytic on the upper half plane. The single poles at ±u!ab — i'J 
correspond to the atomic two-level resonance. On the real frequency axis, x is 
peaked at ±u!ab with the spectral line width 7. Figure D illustrates how the real 
and the imaginary part of the spectral susceptibility depend on the frequency. The 
medium is most dispersive near the resonance frequency Uab where, unfortunately, it 
is also most absorptive. Instead of slowing down hght, the medium turns completely 
opaque. 



4 Electromagnetically-Induced Transparency 

Electromagnetically-Induced Transparency (KIT) has served as a method to 
slow down hght significantly |]I| or, ultimately, to freeze light completely |^ . Like 
other successful techniques, EIT is based on a simple idea fll]]: A control beam of 
laser light couples the upper levels of an atom, and, in this way, the beam strongly 
modifies the optical properties of the atom. In particular, the coupling of the ex- 
cited states affects the transition from the atomic ground state to one of the upper 
states, i.e., the ability of the atom to absorb probe photons with matching transi- 
tion frequency. Destructive quantum interference between the paths of the transition 
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Figure 1: Spectral susceptibility of light in an ordinary dielectric medium. The 
figure shows the real part (solid line) and the imaginary part (dashed line) of x{^) 
as a function of the detuning A = u — ujq in arbitrary units. The function is given 
by equations and (|2^). The line width 7 was set to 1/3 in the units used. 
The imaginary part of x{^) has a peak at the atomic resonance frequency (for zero 
detuning). Outside the resonance the real part grows monotonically corresponding, 
according to equation ([T^, to a positive group velocity. Near the resonance the 
dispersion reaches a maximum. At the resonance the real part of x{^) decreases 
sharply (anomalous dispersion leading to a negative group velocity). However, the 
interesting spectral region of low or negative group velocity is totally overshadowed 
by absorption. 
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Figure 2: Three-level atom in a regime of Electromagnetically-Induced Trans- 
parency. The control beam couples the levels 2 and 3, which inHuences strongly 
the optical properties of the atom for a weaker probe beam tuned to the transition 
1^3. 



process turns out to eliminate absorption at exact resonance [|rT|. A medium com- 
posed of such optically-manipulated atoms is transparent in a spectral region where 
it would otherwise be completely opaque. In the vicinity of the transparency fre- 
quency uq the medium is highly dispersive, i.e., the refractive index changes within 
a narrow frequency interval. In turn, probe light pulses with a carrier frequency at 
c^o travel with a very low group velocity Vg [l^ . 

Consider the three- level atom illustrated in figure 0. The atom is characterized by 
the energy- level differences huu and hLJ23 with + UJ23 = = uq. Typically, the 
transition frequencies cuis and UJ23 are in the optical range of the spectrum or in the 
near infrared (lO^^Hz), whereas the frequency uu is much lower (lO^Hz). The atom 
is subject to fast relaxation mechanisms (lO^Hz) that transport atomic excitations 
from the | 3 ) state down to 1 1 ) and from | 3 ) to | 2 ) , mainly caused by spontaneous 
emission. Hardly any excitations move from | 2 ) to 1 1 ) , because the spontaneous 
emission rate is proportional to the cube of the frequency [jl3[. Here the relaxation 



may be dominated by other processes, for instance by spin-exchanging collisions. 
Without relaxation the dynamics of the atom is governed by the Hamiltonian 



H 













\ '^13 
\ '^23 
Uu)x3 



(30) 



The Hamiltonian represents the atomic level structure and describes the —PE in- 
teraction with light, considering here only the frequency components E^"* and E^^ 
that match approximately the level structure. The Ep and Ec fields are the probe 
and control light respectively. We describe relaxation phenomenologically by the 
flip processes 

ii = |l)(3|, i2 = |2)(3|, (31) 
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occurring at the rates 71 and 72, typically a few 10® Hz. Suppose that the three- 
level atom is illuminated with monochromatic control light at frequency uj^ in exact 
resonance with the 2^3 transition, 



1^23 • 



(32) 



Consider a regime of linear response. In this case we can decompose the probe 
field into monochromatic waves, to describe completely the reaction of the atom. 
We characterized the two light fields involved by their Rabi frequencies VL^ and 
defined as 



Vtc e 



(33) 



The Rabi frequencies set the time scales of atomic transitions caused by the applied 
light fields. The control beam shall dominate all processes. 



^lc\ > I ^^pl , 7i , 72 



(34) 



Mediated by relaxation, the atomic dipoles lose any initial oscillations they might 
have possessed and follow the optical fields. To find the stationary state, we utilize 
an interaction picture generated by 





uj — Uc 
CO 



(35) 



Due to commutation relations of the type ( pOD the dissipative part of the master 
equation is not changed in the interaction picture. The transformed Hamiltonian 
has become time-independent. 



H = -h 











2 "p 



co-ujo i ^]* 



2 "^'p 



(36) 



Similar to a two-level atom in a stationary state (|2l 
balance the relaxation processes. 



the optical transitions should 



^ [p, H] = j:ii {^lAp - 2ArpA\ + pA\A, 



(37) 



We could solve exactly the linear equation (^) for the matrix elements of p with 
trp = 1 (using computerized formula manipulation, for example), but without gain- 
ing much insight. Fortunately, since we are interested in the regime (|3^), we can 
find transparent approximations. Suppose first that also 



(3^ 
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We expand the solution of equation 
Ml) and (|38|), and get 



|37|) to quadratic order in the small quantities 





-. 1 


n; 


p = 




1 pi 














2(^ - u;o) 






"p 



(39) 



We proceed similarly to our analysis of the two-level atom and find, in the positive- 
frequency range, the spectral susceptibility 



2a , 
X= — 



given here in terms of the parameter 



a 



2 S{)\Ec\ 



(40) 



(41) 



The spectral susceptibility x depends linearly on the detuning oj — ojq and vanishes 
at the resonance frequency. Here the phase velocity of light is exactly the speed of 
light in vacuum, c, but the group velocity ([T7|) is reduced by (1 + a), 

(42) 



1 + a 



We call the a parameter ( ^T]) group index. The parameter is proportional to the 
ratio between the probe-hght energy per photon, Uujq, and the control-light energy 
per atom, eQ\Ec\'^ /ua- Consequently, the less intense the control beam is the slower 
the probe light is, a paradoxical behavior. Taken to the extreme, the group velocity 
would vanish when the control beam is totally dimmed. However, in the stationary 
regime that we are considering, the control beam should dominate ([3^ ) and the 
detuning should be small compared with the modulus of the control's Rabi frequency. 

Apparently, the linear spectral slope ( ^OD of the susceptibility ought to be limited. 
To find the limitation, we expand the exact stationary state of the master equation 
(p2D to linear order in f2p, in accordance with a regime of linear response. We obtain 
a spectral susceptibility (^) with 



ri3 



(u - ujq)^ + i{uj - cJo)(7i + 72) - 7 l^c 



(43) 



One can easily verify that the poles of x ^-re located on the lower half plane of the 
complex frequency u;, in agreement with the causality of the medium's response. We 
expand formula ( ^3|) in powers of {uj — ujq) /\VL(.\ and see that the medium becomes 
dissipative when the condition 



7i + 72 



(44) 
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Figure 3: Susceptibility of the probe light in a medium with Electromagnetically- 
Induced Transparency. The figure shows the real part (solid line) and the imaginary 
part (dashed line) of the spectral susceptibility x{uj) as a function of the detuning 
A = uj — ujq in units of the Rabi frequency Qc of the control beam, given by equations 
( |2^ ) and (^^. The line width 71 + 72 was set to r2c/3. The parameters used agree 
with the ones in figure Comparing the two figures, we see that EIT radically 
alters the susceptibility in a spectral region around the probe resonance uq. Here 
the imaginary part of x{uj) vanishes. The medium has become transparent where 
it would be completely opaque without the influence of the control beam. In the 
transparency window the real part of x{^) increases linearly with a steep slope, 
indicating that the medium is extremely dispersive. As a consequence of equation 
( |r^ ) the group velocity of the probe light is significantly reduced. 



is violated. For a large detuning we can ignore the | |ficP term in the susceptibility 
(^Sp. We obtain the simple Lorentzian of an ordinary medium, with the com- 
bined line width 7 = 71 + 72- Outside the narrow transparency window of EIT, the 
absorption of the medium has even slightly increased, because the control beam cou- 
ples the medium atoms to a second dissipative transition process. The maximally 
tolerable detuning for transparency is proportional to the group velocity, since Vg is 
proportional to | f^cP for sufficiently slow light. In practice the detuning is usually 
limited by eVgUJo/c with e in the order of a few 10^^. The transparency window 
concerns slow light in moving media, because of the Doppler effect. An atom with 
velocity u causes a Doppler detuning of uujq/c. If we fix the spectral range in the 
laboratory frame, the maximally tolerable velocity is evg. EIT is velocity-selective. 
Figure |^ illustrates the spectral susceptibility. 
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5 Dark-state dynamics 



Suppose that a dominant and monochromatic control beam has, after relaxation, 
prepared the atom in the stationary state How will the atom evolve when 



the control and probe strengths vary |14|? First we note that the state (pOD is 



statistically pure [to quadratic order in the small quantities (|3^) and (^8|)] so that 

tr{p2} = l. (45) 



When the purity condition (^5|) is satisfied the density matrix contains a single state 
vector 

P=|V^o)(^o| (46) 

with, in our case, 

l*o> = C>.iVo(|l)-||2> + H<^n,|3)). (47) 

The stationary state does not depend on the relaxation rates but only on the pa- 
rameters of the Hamiltonian (^). Remarkably, even when the parameters vary, the 
state is protected from further relaxation, as long as the | 3 ) component is small, 

P33 = (3|p|3)«l. (48) 

Once the atom is in a pure state with sparsely populated top level, the purity (|45|) 
does not change during the evolution (^, 

dtr{/52} = 2tr{pdp} = 4 [71(1 - pn) + 72(1 - P22)] P33 dt . (49) 

The pure state so protected is called a dark state |jl5[. Although dark states are 



initially prepared due to the relaxation of the atomic dipoles, having so adapted to 
the light fields, they are no longer prone to dissipation. 

Suppose that the control and the probe strengths vary. How does a dark state 
follow the light? In the case ( ^Sf ) the state of the atom is dominated by its compo- 
nents in the subspace spanned by the two lower levels | 1 ) and | 2 ). If we find a state 
I ip ) that describes correctly the dynamics (^ in this subspace, the third component 
(3|'0) must be correct as well, to leading order in P33. The lower ranks enslave 



the top level. Since the relaxation processes (|3TD do not operate within the lower 
subspace, we can ignore dissipation entirely, to find the dominant state of the atom. 
We represent both control and probe light in terms of variable Rabi frequencies, 

g-.-.* ^ ^ , np e-'^^' = ^ , (50) 

defined here with respect to the atomic transition frequencies = UJ23 and ujq = uji^. 
We write down the state vector 

\,^^U.,nL)-^12)^'M,8.%3)) (51) 
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with the abbreviations 



Uo 
N 



1 

e 






iuji2t 






-iujQt 



je 



1 + 



-1/2 



In a stationary regime under the conditions ( plf ) and 
with the dark state (P7|). We see from the properties 



(52) 
(53) 
(54) 
(55) 

the vector (|5lD agrees 



(56) 



that I ip ) satisfies the differential equation 

ihdtlip) = H\iIj) +thdt {31-^)13) . 



(57) 

Consequently, the vector (^T]) describes correctly the dynamics of the atom in the 
lower-level subspace. Therefore, the atom remains in the dark state (0), as long 
as the atom's evolution never leads to an overpopulation at the top level | 3 ). The 
initial relaxation-dominated regime has prepared the dark state, but later the atom 
follows dynamically without relaxation ||14|| . 

We calculate the matter polarization ([T0|) generated by the dark states of the 
atoms that constitute the medium. The positive-frequency component of P is 



p(+) 



^K3l(3|^)(^|l) 



2 



K31 e 



K 



31 







n lap 



idt -ujQ-i 



idt\n,\) 



+ 0C 



(5^ 



with 6c = argfic- Assume, for simplicity, that Qc is real. Otherwise we can easily 
incorporate the phase 6c of the control field in the phase of the electric field without 
affecting the wave equation (^), as long as 6c varies slowly compared with the optical 
frequency Uq. We adopt the definition ( Pl) of the group index a, and get 



£0 



. a 

''2a 



p 



We approximate 
2u}o{idt - 



(59) 



(60) 
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and obtain from the general wave equation (^) an equation that is vahd for both 
the positive and the negative frequency component of the probe hght, 



dl-^dl + N^{dtadt + 



E„ = 0. 



(61) 



The dark-state dynamics may lead to a non-linear effect of the medium, described 
by the factor in the wave equation (|6ll). However, when the probe is significantly 
weaker than the control light, the medium responds linearly, 

(62) 



dt{l + a)dt - c^di + auj^jEp = 0. 

This wave equation governs the propagation of slow light in a regime of linear re- 
sponse and undisturbed dark-state dynamics. 



6 Slow-light pulses 

Consider a pulse of probe light in an EIT medium with variable group index (^if ). 
Suppose that the group velocity (^2]) does not vary much over the scale of an optical 
wavelength (0.5 x 10~®m) or an optical cycle (10~^^s). In this case we could apply 
the Hamiltonian theory (|T6|) of a fictitious light particle to predict the position of 
the pulse peak. Because particle trajectories must not split, a slowly varying group 
index cannot cause reflection. Suppose that the pulse is traveling to the right. Then 
the pulse will continue to do so, and we can express the slow-light wave as 

Ep = £ expHkz — iut) + S* exp(—ikz + iut) , k = — , (63) 

c 

assuming that the electric-field amplitude S is slowly varying compared with the 
rapid optical oscillations. We approximate 



exp{-ikz + iut) d'^E^^^ ^ [-u^ - 2iujdt) £ , 
exp(— ■jA;^; + iut) dtE^^ ^ —iuj8 , 

ex^{-ikz + iojt) dlE'^^ ^ (-k'^ + 2ikdt) 8 . (64) 



and get from the wave equation ([6^ ) 

-2iuo{l + a)dtS ^ ({l + a)uo'^ + iua-c^k'^ + 2ikc^d,-auo'^S 



= (2iujcdz + iwa + 0(00'^ — cUq)) S 

= 2icu (^cdz£ + I ^) (65) 

when the carrier frequency 00 is equal to the transparency resonance Uq. We thus 
obtain the propagation equation |]T^ 

(1 + a) dtS + cdzS + ^£ = (66) 
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In order to understand the principal behavior of ordinary slow-hght pulses, we con- 
sider two cases — a spatially varying yet time-independent group index and a spa- 
tially uniform yet time-dependent a. 

When the group index does not change in time, the propagation equation (|66|) 
has the simple solution 

E{t,z) = E,{t- /-), (67) 

J Vg 

in terms of the group velocity (^21). At each space point z the pulse raises and falls 
in precisely the same way. However, because light is slowed down, the spatial shape 
of the pulse shrinks by a factor of Vg/c compared with the pulse length in vacuum, 
for example by 10"'' for a group velocity of 30m/s. The intensity of the pulse is 
unaffected, despite the enormous pulse shortening, and the pulse energy has gone 
into the amplification of the control beam. 

Consider the other extreme, a spatially uniform EIT medium with adjustable 
group velocity. In this case, the solution of the propagation equation (|66D is 



£{t,z) = £o(^z - J VgdtJ \lvg/c. (6^ 



The pulse envelope £ propagates again with the group velocity Vg but the pulse 
length is not changed. However, the spectrum of the pulse around the carrier fre- 
quency ujq shrinks by a factor of Vg/c. Additionally, the intensity drops by Vg/c as 
well. The ratio between the control (5l) and the pulse intensity (|68D remains large. 



even in the limit of a vanishing control field when Vg vanishes, as long as is 
large (for a sufficiently dense medium). Therefore, the reaction of the probe field is 
remarkably consistent with the requirements for undisturbed dark-state dynamics 
[|T^. One can freeze light without losing control 0, |^. 

7 Effective field theory 

After having studied two examples of pulse propagation in an EIT medium, we 

develop an effective field theory of slow light. We generalize the wave equation ( |H^ ) 
to three-dimensional space and calibrate the electric field strength in appropriate 
units, 

Ep(t,x)=(-j ^(t,x). (70) 
We introduce the Lagrangian 

^ = ^ ((1 + amvf - c'iV^r - ao^oV) (71) 

and see that the wave equation (|62|) is the resulting Euler-Lagrange equation. We 
have chosen the prefactor of the Lagrangian (|7lD such that ^ agrees with the free- 
field Lagrangian (|I|) for zero a and frequencies around c^o- Therefore we regard Jif 
as the effective Lagrangian of slow light. 
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Let us use the Lagrangian (^) to calculate the energy balance of slow light. 
According to Noether's theorem |]16| we obtain the energy density 

/ = !((! + aM^r + c'iV^r + a cUoV') (72) 
and the energy flux (Poynting vector) 

S = -hc\dt^){V^). (73) 
As a consequence of the wave equation (|62D we obtain the energy balance 

dJ + V-S = ^{^' + u;^) . (74) 

Temporal changes in the control field, modifying the group index (^), do not con- 
serve energy. In fact, the experiment |0] indicates that the control beam can amplify 
light stored in an EIT medium with zero group velocity. In the experiment , slow 
light enters the EIT sample and is then frozen inside by turning off the control field. 
Switching on the control releases the stored light. The pulse emerges with an inten- 
sity that depends on the control field and that may exceed the initial intensity, in 
agreement with equation ([681) . Clearly, this phenomenon is only possible if energy 
is indeed transferred from the control beam to the probe light. 



8 Moving media 

An EIT medium slows down light, because the medium is extremely dispersive, re- 
acting differently to the different frequency components of a pulse. The peak position 
of the pulse is the place where the frequency components interfere constructively. 
By slightly modifying the phase velocity of each component the medium influences 
strongly their interference, slowing down the pulse dramatically. 

The extreme spectral sensitivity of slow light can be also applied to observe 
optical phenomena in moving media, caused by the Doppler effect. A uniformly 
moving medium would not present an interesting case, though, because uniform 
motion just produces a global frequency shift. However, slow light is a superb tool 
in detecting non-uniform motion such as rotation ||T^. To understand the principal 



effect of slow light in moving media, we modify the Lagrangian (^T]) to account 
for the Doppler effect. We assume that (^) is valid in frames co-moving with the 
medium and transform back to the laboratory frame. We note that {dt^pY — (V ip)"^ 
is a Lorentz invariant and focus on the first-order Doppler effect in the a{dtipy term, 
assuming the realistic case of non-relativistic medium velocities. We replace dt by 
dt + u-V and neglect the term quadratic in u. In this way we obtain the Lagrangian 
of slow light in a moving medium 

^ = ^o-^S-u, S = -hc\dty^)iV^) (75) 

in terms of the Lagrangian for the medium at rest (0). We see that the flow u 
couples to the Poynting vector S of slow light, similar to the Rontgen interaction of 
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moving dipoles in electromagnetic fields |T^. We obtain from the Euler-Lagrange 
equation 



(dt{\ + a)dt - c^V^ + a cUq + 9ta u ■ V + V ■ au 9; 







(76) 



with the differential operators acting on everything to the right. For frequencies 
near the EIT resonance o^o we represent the positive-frequency part <y9'^+) of ip as 



(77) 



and perform similar approximations as in Section 5. We obtain the Schrodinger-type 
equation 



A 



1 



2v„ 



a 



V^o = t: -^^V + - u ifQ 



^9 \ ^"^9/ 

with the effective Planck constant reduced by 2ti 

A = ^. 



(7J 



(79) 



The flow has a two-fold effect: On one hand, the velocity u appears similar to an 
effective vector potential, for example as the magnetic vector potential acting on 
an electron wave, and, on the other hand, the hydrodynamic pressure proportional 
to acts similarly to an electric potential. A vortex flow will generate the optical 
equivalent of the Aharonov-Bohm effect, see reference [|l^] for details. 

A moving slow-light medium is also equivalent to an effective gravitational field 
[|]. Consider monochromatic light at exact resonance frequency cuq- Iii this case, 
we can write the wave equation (|76D in the form of a Klein-Gordon equation in a 
curved space-time. 



.if 



{+) 



0, 



with = {dt/c,V) and 



fJ,U 



-9 9' 



flU 



1 au/c 
au/c —1 



^0) 



(81) 



Here g'^'^ represents the effective space-time metric experienced by monochromatic 
slow light in a moving medium, to first order in m/c. We easily find the determinant 
g of the inverse of g'^'^ from the relation 



det/ = -g^g = -g = - {l + a'- 



u 



The effective space-time line element ds^ is, up to a conformal factor, 

ds^ = dt^ + 2a dt u • dx - dx^ , 
resembling the line element of a moving coordinate system. 



(82) 



^3) 



(c^ - u^)dt^ + 2 dt u ■ dx - dx^ 



^4) 
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for example of a rotating system. The parameter a quantifies the degree to which 
the motion of the medium is transferred to the propagation of hght in the medium, 
the degree of dragging. The group index is thus equivalent to Fresnel's dragging 



coefficient ||T^ For slow light a is very large indeed. Therefore, slow light is 
able to sense minute flow variations. Even subtle quantum flows imprint phase shifts 
onto slow light that are detectable using phase-contrast microscopy . 



Sound waves in a fluid experience the flow as an effective space-time metric as 
well [^. The acoustical line element is proportional to the element (|^) of a moving 



system of coordinates, with two crucial differences: The flow is not subject to the 
rigidity of moving coordinate systems and, more importantly, in the acoustic metric 
the speed of light, c, is replaced by the speed of sound. A supersonic flow surpasses 
the sound barrier and can, under suitable circumstances, generate an artificial event 
horizon where the fiow speed u reaches c. Here it is necessary that the goo element 
of the metric vanishes. In the slow-light metric the all-important term —u^dt"^ 
is missing in ^fQO, at least to the level of approximation we are considering here. 
We obtain a term proportional to —v?(\.t^ when we include effects of higher-order 
Doppler detuning. The critical velocity is then of the order of c/ ^/a. To observe the 
quantum effects of light generated by a horizon we would need to employ a steep 
profile of the fiow speed. This causes a severe problem, because the Doppler detuning 
will exceed the transparency window of EIT. The Doppler effect plays a beneficial 
role in the sensitivity of slow light to motion, but it will also cause significant light 
absorption when one attempts to reach an artificial event horizon. The medium 
will certainly turn black, but not into a black hole. However, one could employ a 
spatially varying profile of the group index to create an interface that resembles an 
event horizon for slow light and that avoids this problem pB[. Slow light offers a 



variety of options for interesting experiments exploiting the analogs of light in media 
with effects in other areas of physics, and new ideas continue to emerge. 



9 Summary 

Light has been slowed down dramatically ||T| or even stopped completely 0, To 
understand how this feat has been achieved, we studied the physics behind slow 
light, starting from basic first principles of the light-matter interaction. We first 
turned to ordinary optical media, so as to later contrast them with slow-light media 
based on Electromagnetically-Induced Transparency. We studied slow light in two 
regimes — in a stationary regime both dominated and limited by relaxation and 
the control light, and in a dynamic regime almost free from dissipation. Then we 
analyzed the typical behavior of slow-light pulses, before developing an effective field 
theory of slow light that we have subsequently generalized to moving media. 

It is certainly amazing how much a clever combination of atomic physics and 
optics can achieve, but it is also important to understand the principal limits of the 
techniques applied. These limits depend on the details of the physics behind the 
scene. We have tried to elucidate the details of slow light without going into too 
much detail, using models that are simple, but not too simple. 
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